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On the Origin of Strong-Field Polarity Inversion Lines 



B. T. Welschi and Y. Li^ 

Abstract. Several studies have correlated observations of impulsive solar 
activity — flares and coronal mass ejections (CMEs) — with the amount of 
magnetic flux near strong-field polarity inversion lines (PILs) in active regions' 
photospheric magnetic fields, as measured in line-of-sight (LOS) magnetograms. 
Practically, this empirical correlation holds promise as a space weather forecast- 
ing tool. Scientifically, however, the mechanisms that generate strong gradients 
in photospheric magnetic fields remain unknown. Hypotheses include: the (f) 
emergence of highly twisted or kinked flux ropes, which possess strong, opposite- 
polarity fields in close proximity; (2) emergence of new flux in close proximity 
to old flux; and (3) flux cancellation driven by photospheric flows acting fields 
that have already emerged. If such concentrations of flux near strong gradients 
are formed by emergence, then increases in unsigned flux near strong gradients 
should be correlated with increases in total unsigned magnetic flux — a sig- 
nature of emergence. Here, we analyze time series of MDI line-of-sight (LOS) 
magnetograms from several dozen active regions, and conclude that increases 
in unsigned flux near strong gradients tend to occur during emergence, though 
strong gradients can arise without flux emergence. We acknowledge support 
from NSF-ATM 04-51438. 



1. Strong gradients across PILs 

It has been known for decades that flares and filament eruptions (which form 
CMEs) originate along polarity inversion lines (PILs) of the radial photospheric 
magnetic field. In studies using photospheric vector magnetograms. Falconer 
et al. (2003, 2006) reported a strong correlation between active region CME 
productivity and the total length of PILs with strong potential transverse fields 
(> 150 G) and strong gradients in the LOS field (greater than 50 G Mm~^). 
They used a ib2-day temporal window for correlating magnetogram properties 
with CMEs. Falconer et al. (2003) noted that these correlations remained essen- 
tially unchanged for "strong gradient" thresholds from 25 to 100 G Mm~^. Using 
more than 2500 MDI (LOS) magnetograms, Schrijver (2007) found a strong cor- 
relation between major (X- and M-class) flares and the total unsigned magnetic 
flux near (within ~ 15 Mm) strong-field PILs — defined, in his work, as regions 
where oppositely signed LOS fields that exceed 150 G lie closer to each other 
than the instrument's ~ 2.9 Mm resolution. Schrijver's (2007) effective gradient 
threshold, ~ 100 G Mm~^, is stronger than that used by Falconer et al. (2003, 
2006). 

Although these studies were published recently, the association between 
flares and 5 sunspots, which posses opposite-sign umbrae within the same penum- 
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bra — and therefore also possess strong-field PILs — has been well known for 
some time (|Kiinzel 1960|;_l s ammi s et al. 2000l '). In particular, [3^5 spot groups 
are most likely to flare ( Sammis et al. 20001 ) . A 0'y designatio n means no obvi- 
ous north-south PIL is present in an active region (|Zirin 19881 ). 

We note that Cui et al. (2006) found that the occurrence of flares is corre- 
lated with the maximum magnitude of the horizontal gradient in active region 
LOS magnetograms — not just near PILs — and that the correlation increases 
strongly for gradients stronger than ~ 400 G Mm~^. 

One would expect the measures of CME- and flare- productivity developed 
by both Falconer et al. (2003,2006) Schrijver (2007) to be larger for larger ac- 
tive regions. Importantly, however, both studies showed that their measures of 
flux near strong-field PILs is a better predictor of flare productivity than total 
unsigned magnetic flux. Evidently, more flux is not, by itself, as significant a 
predictor of flares as more flux near strong-field PILs. 

These intriguing results naturally raise the question, "How do strong- field 
PILs form?" For brevity, we hereafter refer to strong-field PILs as SPILs. 

Schrijver (2007) contends that large SPILs form primarily, if not solely, by 
emergence. But he also noted that flux emergence, by itself, does not necessarily 
lead to the formation of SPILs. Rather, a particular type of magnetic structure 
must emerge, one containing a long SPIL at its emergence. He suggests such 
structures are horizontally oriented, filamentary currents. 

Beyond the "intact emergence" scenario presented by Schrijver (2007), other 
mechanisms can generate SPILs. When new flux emerges i n close proximity to 
old flux — a common occurrence (iHarvev and Zwaan 19931 ) — SPILs can form 
along the boundaries between old and new flux systems. Converging motions in 
flux that has already emerged can also generate SPILs. If the convergence leads 
to flux cancellation by some mechanism — eme rgence of U lo ops, submergence 
of inverse-U loops, or reconnective cancellation ( Welsch 20061 ) — then the total 
unsigned flux in the neighborhood of the SPIL might decrease as the SPIL forms. 
We note that, while cancellation in already-emerged fields can occur via fiux 
emergence (from upward moving U-loops), the emergence of a new flux system 
across the photosphere must increase the total unsigned flux that threads the 
photosphere. 

If the emergence of new flux were primarily responsible for SPILs, then 
a straightforward prediction would be that an increase in total unsigned flux 
should be correlated with an increase in the amount of unsigned flux near SPILs. 
Hence, observations showing that increases in the unsigned flux near SPILs 
frequently occur without a corresponding increase in total unsigned flux would 
rule out new flux emergence as the sole cause of these strong field gradients. 

Our goal is to investigate the relationship between increases in the amount 
of unsigned flux near SPILs with changes in unsigned flux in the active regions 
containing the SPILs, to determine, if possible, which processes generate SPILs. 



2. Data 

From days-long time series of deprojected, 96-minute, full-disk MDI magne- 
tograms for Nar = 64 active regions, we computed the rates of change of un- 
signed flux near SPILs, following the method described by Schrijver (2007). We 
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also computed the rates of change of total unsigned line-of-sight magnetic flux 
these active regions. 

Our active region sample was chosen for use in a separate study of the rela- 
tionships between surface flows derived from magnetograms and CMEs. For the 
purposes of that study, we typically selected active regions with a single, well 
defined PIL, for ease in identifying the presence of shear i ng and/or converg - 
ing flows some CME models employ ( Antiochos et al. 19991 : Linker et al. 20dll ). 
The sample used here includes regions from 1996 - 1998, and includes regions 
that did and did not produce CMEs. Some of our magnetograms image the 
same active region as it rotated back onto the disk one or more times. Also, 
some of our selected regions are so decayed that they lack spots, and therefore 
have no NOAA designation. 

Here, we analyze A^mag = 4062 magnetograms. Pixels more that 45° from 
disk center were ignored. To convert the LOS field, -BloS) to an estimated radial 
field, Bfi, cosine corrections were used, Br = Bi^os/ cos(0), where O is the angle 
from disk center. 

Triangulation was used to interpolate the Br data — regularly gridded in 
the plane-of-sky, but irregularly gridded in spherical coordinates {6, cj)) on the 
solar surface — onto points {9', (j)') corresponding to a regularly gridded, Merca- 
tor projection of the spherical surface. This projection was adopted because it 
is conformal (locally shape-preserving), necessary to ensure displacements mea- 
sured in the tracking study mentioned above were not biased in direction. For 
computing gradients, a conformal projection is also appropriate. The back- 
ground grayscale in Figure [1] is a typical reprojected magnetogram. We note 
that the price of preserving shapes in the deprojection is distortion of scales; 
but this can be easily corrected. 

Each active region was tracked over 3-5 days, and cropped with a moving 
window. A list of tracked active regions and mpeg movies of the active regions, 
are online, at http://sprg.ssl.berkeley.edu/~yanli/lct/ . 



3. Analysis methods 

To identify SPILs, we used the gradient identification technique of Schrijver 
(2007). For a magnetogram at time tj, binary positive/negative strong- field 
masks — where Bn > 150 G and Bji < —150 G, respectively — were con- 
structed, then dilated by a (3x3) kernel to create dilated positive and negative 
bitmaps, M±. These procedures are illustrated in Figured) Regions of overlap, 
where Mql = M+M_ = 1, were identified as SPILs. In Figured! Mql / for 
a single pixel, at (x,y) = (112,91). 

To quantitatively define neighborhoods around SPILs, Mql is convolved 
with a normalized Gaussian, 

G(u,v) = Go^exp(-[u2 + 7;2]/2a2) (1) 

where Gq = J du J dv ex.p{—[u'^ + f^]/2(T^), and a = 9 pixels (corresponding to 
a FWHM ~ 15 Mm at disk center), to create "weighting maps," Cmg, where 



Cmg{x, y) = convol(MoL(2;, y),G) 



(2) 
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Figure 1. The background grayscale is typical, reprojected magnetogram; 
white is positive flux, black is negative flux. The inner black and white con- 
tours enclose signed, strong-field masks — regions where Bji > 150 G and 
Bn < —150 G (respectively). The outer black and white contours show the 
outlines of M± , dilated bitmaps of the strong- field masks with a 3 x 3 kernel 
function. For this magnetogram, the dilated bitmaps overlap at a single pixel, 
at ix,y) = (112,91). 



Figure [2] shows the product of Bpi with such a weighting map. Following Schri- 
jver (2007), we totaled the unsigned magnetic field over the weighting map to 
determine R, a measure of the unsigned flux near SPILs, 

R = J2\Br\Cmg . (3) 

From a sample of more than 2500 MDI magnetograms, Schrijver (2007) showed 
that R is correlated with major (X- and M-class) flares. 
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Figure 2. Br multiplied by a weighting map, Cmg- R is the total unsigned 
field the window, ^ \Br \ Cmg- 

For each of the =1621 magnetograms with R ^ 0, we summed the 
weighted absolute magnetic field in the previous magnetogram, BR(ti-i), using 
the weighting map from tj, to compute the backwards-difference AR, 

AR = Y,{\BR{ti)\ - \BR{ti.i)\) Cmg • (4) 

We also computed the change in summed, unsigned field, 

AB = Y,{\BR{U)\-\Bn{U-i)\) , (5) 

to determine if new flux is emerging or if flux is canceling. If new flux is emerging, 
we expect AB > 0. If flux is canceling, we expect AB < 0. Like Schrijver (2007), 
we have opted to keep R in units of flux density; for simplicity, we also keep AB 
in these same units. 
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When the overlap map Mql for Bfi{ti) is identically zero, R is also zero, 
and AR and AB are not computed. 

In §L] we discussed processes that can cause changes in R. What processes 
can lead to AB ^ 0? Emergence of new flux or cancellation (both only happen 
at PILs) can make AB ^ 0, and these processes are probably related to evolution 
in R. Flux can also cross into or out of the cropping window. Since our cropping 
windows were selected to include essentially all of each tracked active region's 
flux, systematic errors arising in this way are expected to be small. A more 
severe effect is the "unipolar appearance" phenomenon characterized by Lamb 
et al. (2007), who found that the majority of newly detected flux in the quiet 
sun is due to coalescence of previously existing, but unresolved, single-polarity 
flux into concentrations large and strong enough to detect. While it is unclear 
if the conclusion reached by Lamb et al. (2007) for the quiet sun also applies 
in active regions, this is plausible. Moreover, much as flux can "appear," flux 
can also disappear, via dispersive photospheric flows or perhaps even molecular 
diffusivity. Also, simultaneous emergence of new flux and cancellation of exist- 
ing flux can occur within the same active region, masking the effects of both 
processes. Practically, therefore, we can only refer to increases in unsigned flux 
as "possible new flux emergence," and to decreases in unsigned flux as "possible 
cancellation." 

4. Results and conclusions 

In Figure O we show a scatter plot of changes in i? as a function of changes in 
B. The plot does not show the full range in AB, but the AR for outliers on the 
horizontal axes are near zero. One striking feature of the plot is its flatness, i.e., 
that most changes in B are not associated with any change in R. In Table [H 
we tabulated the data points in each quadrant of this plot. Clearly, increases in 
R, the unsigned flux near SPILs, usually occur simultaneously with increases in 
the unsigned flux over the entire active region. Increases in R only occur less 
frequently when flux is decreasing, i.e., during cancellation. 

Table 1. Breakdown of Flux Changes 
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We set out to answer the question, "How do strong-fleld PILs form?" We 
related changes in total, unsigned flux over whole active regions with changes 
in total, unsigned flux in subwindows of the same active regions — defined by 
weighting maps. One might expect, therefore, that these quantities should be 
correlated, casting doubt about our ability to discrimintate between changes 
in total flux in active regions and in subwindows. If the two were strongly 
correlated, the excess of events with AR > and AB > might not be very 
meaningful. In fact, however, AR and AB are poorly correlated: the two have a 
linear correlation coefficient r = 0.29, and a rank-order coefficient of 0.36. This 
suggests that the relationship between increases in R and increases in total, 
unsigned active region flux is not an artifact of our approach. 
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Figure 3. A scatter plot of changes in i? as a function of changes in B. 
Increases in R, the unsigned flux near SPILs, usually occur simultaneously 
with increases in the unsigned flux over the entire active region. Increases in 
R only occur rarely when flux is decreasing, i.e., during cancellation. For a 
breakdown of the data points in each quadrant, see Table [TJ 



Nonetheless, our active region sample is not ideally suited to address the 
origin of SPILs, generally. Our sample was not unbiased with respect to active 
region morphology; we selected regions with well-defined PILs. In addition, 
our sample included some decayed active regions that NOAA AR designations. 
Consequently, we believe that a follow-up study, with a much larger, unbiased 
sample of active regions, is warranted. 

With caveats, therefore, our study supports Schrijver's (2007) contention 
that the emergence of new flux creates the strong-field polarity inversion lines 
that he found to be correlated with flares. 

Acknowledgments. We acknowledge the support of NSF Grant NSF-ATM 
04-51438. 
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